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ABSTRACT 

TeV gamma rays emitted by GRBs are converted into electron-positron pairs 
via interactions with the extragalactic infrared radiation fields. In turn the pairs 
produced, whose trajectories are randomized by magnetic fields, will inverse 
Compton scatter off the cosmic microwave background photons. The beamed 
TeV gamma ray flux from GRBs is thus transformed into a GeV isotropic gamma 
ray flux, which contributes to the total extragalactic gamma-ray background 
emission. Assuming a model for the extragalactic radiation fields, for the GRB 
redshift distribution and for the GRB luminosity function, we evaluate the contri- 
bution of the GRB prompt and scattered emissions to the measured extragalactic 
gamma-ray fiux. To estimate this contribution we optimistically require that the 
energy fiux at TeV energies is about 10 times stronger than the energy fiux at 
MeV energies. The resulting gamma-ray diffuse background is only a small frac- 
tion of what is observed, allowing blazars and other sources to give the dominant 
contribution. 

Subject headings: 98.70.Rz 



Introduction 



The nature of the extragalactic gamma ray background emission has been a topic of 
great inte rest since EGRET colla boration evaluated its spectrum in the range from 30 MeV to 
100 GeV (ISreekumar et al.lll998l ). The diffuse emission coming from beyond the galaxy was 
determined by subtracting the contributions of resolved point sources, the diffuse galactic 
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emission, and the instrumental background from the gamma-ray intensities observed by 
EGRET. The emission is found to be a power law in energy 



dA dt dn dE 



(7.32 ±0.34) X 10' 



E 



-2.10±0.03 



•0.451GeV^ 



cm^^s^^sr^^GeV"^ 



and is highly isotropic on the sky. iMukherjee fc Chiang] (119991 ) suggested that blazars can ex- 
plain up to 25 per cent of the extragalactic emission. The hypothesis that the flux is predomi- 
nantly due to blaza rs seems to be reinforced by a new evaluation of the extr agalactic emission 



( Strong et al.l 120041 ) which is slightly lower and steeper than that found by lSreekumar et al. 
( 119981 ). The result is not consistent with a power law an d shows some positive curva ture 



which the authors relate to an origin in blazars emission. iKneiske fc Mannheim! (120051 ) re- 
cently sugges ted that up to 85 per cent of the e xtragalactic emission could arise from blazars. 
According to IStecker fc Salamon I (119961 . l200ll ) blazars can account for the entire extragalac- 
tic 7-ray background observed by EGRET. Above 100 MeV normal galaxies contribute from 
3 to 10 per cent of the obse rved diffuse extragalactic flux. However from the analysis of 
( lErlykin fc Wolfendale Ill995l ) the spe ctrum of norma l galax ies seems to differ from the ex- 
tragalactic diffuse emission spectrum. iDar fc Shaviv I (119951 ) have suggested that the extra- 
galactic diffuse arises from cosmic rays interacting with intergalactic ga s. This contribution 



seems t o disagree with th e diffu se gamma ray spectrum according to IStecker &: Salamon 



( 119961 ). iLoeb fc WaxmaxI (120001 ) suggested fossil radiation from shock ac c elerat ed cosmic 



rays during structure formation as possible contribution. IStawarz et al. I (120061 ) have es- 
timated that inverse Compton scattering of starlight photon fields by the ultrarelativistic 
electrons in kiloparsec-scale jets in PR I radio galaxie s con tributes about one percent to the 
EGRET extragalactic flux. Ichi fc Wolfendale! J 19891 ) and Iwdowczvk fc Wolfendalel Jl99oh 
have shown that upscattering of CMB to 7-ray energies by cosmic ray electrons and protons 
does not provide a sinificant contribution to the diffuse extragalactic emission. 

However, the question about the origin of the extragalactic emission is still open and can 
possibly be solved by admitting that different sources contribute to it. In fact all unresolved 
discrete sources outside the Galaxy contribute to the extragalactic background emission; 
the problem is clearly to unders t and h ow big the different contributions are. As already 
pointed out by iHartmann et al.l (120031 ) . who considered GRBs as a source for the diffuse 
gamma-ray emission at MeV energies, prompt and delayed emissions from GRBs should also 
contribute to the diffuse extragalactic emission, especially if some GRBs emit photons in the 
GeV-TeV energy. In fact, outside the GRB source, due to interactions with cosmic infra-red 
background photons, most of the high energy GRB photons produce high-energy electron- 
positron pairs. The pairs inverse Compton scatter off CMB photons and produce secondary 
photons, which in turn interact with IR photons and generate other pairs. Multiple inverse 
Compton scatterings occur until the energy of the secondary photons is no longer large 
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enough to trigger pair production with the IR photons. When the energy of the scattered 
photons is not sufficient to produce a subsequent pair, the photons travel to the Earth 
without undergoing further absorption. After multiple pair-production and inverse Compton 
processes, the initial e nergy of the TeV photons e scaping the GRBs will have been shifted 



to M eV-GeV energies flCoDPi k Aharonianl ll997l: IPlaeal Il995l : IPai fc Lul 12002 : IWang et al. 



2004 : iRazzaque et al.ll2004j ). The data on extragalactic diffuse emission in the MeV-GeV 
energy range are thus useful to put constraints on high energy emission from gamma ray 
bursts at GeV-TeV energy. 



In the following we assume a flat cosmological model with normalized Hubble constant 
0.71, with V Lrn. = 0.3 and = 0.7. Cosmological distances and volumes are calculated 



h 



following iHoggl (119991 ) 



GRB model 



Gamma ray bursts are cosmological spots of short, intense and narrowly beamed 7- 
ray emission, whose observed isotropic luminosities, Lj^o, are in the range of 10^^ — 10^^ 
erg/s. The huge energy released by GRBs is almost certainly produced by ultra-relativistic 
flows, whose bulk Lorentz factors may vary in the range of 100 and 1000, w ith a typical 
value of 300 jBaring fc Harding|ll997l : Ihithwick fc SarilboOll : Izhang et al.lbooeh . The widely 
accepted ffieball internal-external shock model for GRBs envisages that high energy protons 
and electrons are accelerated in shocks by interaction with magnetic fields. The process of 
Fermi acceleration leads to electron and proton power-law spectra of index between -2 and 
-3. Electrons and protons cool through synchrotron emission. Since the electron cooling 
time is shorter than the proton cooling time, as a ffist approximation, the GRB prompt 
spectrum is believed to be dominated by electron synchrotron emission only. T he observed 



prompt spectrum arising from electron synchrotron emission is a Band function (IBand et al. 
19931 ). proportional to E°' below the synchrotron peak energy, E^k, with a about -1 and to 



above the peak energy, with j3 usually between -2 and -3. Below the synchrotron self 
absorption energy, Egsa-, low energy photons are absorbed by ffieball electrons in a magnetic 
field by synchrotron self-absorption. 



^qO ^q2 

Energy (GeV) 



Fig. 1. — Band spectrum and Inverse Compton spectrum for a GRB having isotropic lumi- 
nosity 10^^ erg/s and located aX z = 0.1. The slope a = —1 and (3 varies bewteen -2 and 
-3 for both synchrotron and IC spectra. The ratio r between the synchrotron peak flux and 
the TeV emission peak flux is assumed to be 10. 
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High energy (up to TeV) emission has been detected in some GRBs ( IGonzales et al. 
2OO3I : iHurlev et aP Il994j : ktkins et all boosh. According to the fireball mod el, TeV pho- 



tons can be produced in both internal (jPai fc Lul 120021 : iRazzaque et al.l |2004J ) and external 
( Zhang &: Meszaros 2001 : Bottcher fc Dermerl IQQsi) shocks, either through electron inverse 
Compton or proton synchrotron emission. IWang et al.l (120011 ) proposed that the cross inverse 
Compton scattering between the photons and electrons in forward and reverse shocks can 
also produce TeV photons. Tev photons can be also generated through inverse Compton 
scattering of shock accelerated particles in external shocks off a bath of photon that overlaps 
the shocked region. The photon bat h could be either the prompt gamma-ray emission itse lf 
( iBeloborodovl l2005l : iFan et al.ll2005l ) or the late-time X-ray flares ( IWang fc Meszarosll2006l ). 
For such distant sources as GRBs, TeV photons are mostly absorbed through pair production 
with cosmic background radiation (GBR), but some indications of TeV gamma ra ys from 
GRBs at low redshi ft were provi ded by experiments like M ilagro (lAtkins et al.ll2000l ). Hegra 
( jPadilla et al.l Il998l ) and Tibet (lAmenomori et al.l 119961 ). In our calculation TeV prompt 
energy photons are taken into account in the prompt spectrum without investigating the 
details of the way they are created. We model the GRB prompt TeV emission as an addi- 
tional broken power law with low energy spectral index a = —1 and high energy spectral 
index = —2. The syn chrotron peak energy Epk depends on the isotropic energy of the 
burst (lAmati et al.ll2002l ). In order to match the typical observed Epk at 100 keV - 1 MeV, 
the typical random electron Lorentz factor 7e is about 300-1000. Since most high energy 
emission components are related to inverse Compton (IC) scattering, we design our high 
energy component to mimic the IC in the internal shocks. The typical separation between 
the synchrotron typical frequency and the IC typical frequency is 7^, which varies between 
10^ and 10^. So the break energy of the additional TeV compo nent is chosen to be around 
10^ times the synchrotron peak Epk ( IZhang fc Meszarosll2002al ). The relative energies con- 
tained in the synchrotron component and the IC component depend on t he IC y-parameter 
(IPanaitescu fc Kumai][200ol : ISari fc EsinlboOll : Izhang fc MeszaroslboOll ) 



Y 



'IC 



■'syn 



^1 



(2) 



for ee ^ tBi where ee and es are shock energy equipartition parameters for electrons and mag- 
netic fields, respectively. Because of a possible Klein-Nishina limitation, this Y parameter 
could not be too high. An optimistic value would be around 10, corrsponding to ~ lOOe^ 
(see Fig. [1]). This requires that the GRB internal shocks are not very magnetizeqj. 



Some evidence howev er suggests a magnetized central engine |Zhang et al.l booi IPan et al.l booi 
Kumar fc Panaitescul 2003 ). In that case, the IC component is at most comparable to the synchrotron 



component. 
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GRB emission is likely narrowly beamed. The jet may be top-hat shaped w ith uni- 
form energy distribut ion inside the jet cone and a sharp drop off at the edge (IRhoads 



1999 



Sari et al. 



1991 or " structured" with angle- dependent energy density inside the jet 



Zhang fc MeszarosI l2002bl : IRossI et alibood l The total amount 
rected for the beaming effect is likely standard, about 1.3 x 10^^ erg (iFrail et a 



emission energ y cor- 



average beaming factor ranges in the literature from around 75 (iGuetta et al 
krail et al.lboOlh . 



200 ih . The 



20051) to 500 



Photon-photon absorption through pair production inside the GRB source region con- 
trols the ran ge in energy and the a mount of radiation emitted. Our following treatments 
closely follow iRazzaque et al.l (120041 ). In Fig. [2] we plot the optical depths corresponding to 
the different processes taking place inside the GRB source region, 77 absorption, electron 
Compton scattering and 67 ^ e^. The GRB isotropic luminosity is Liso = lO^^erg/s and 
the observed time variability is 6t = 0.01 s. The time variability is an important parameter 
for the description of gamma ray bursts. In fact the shock radius 



rsh = 2c5tTl 



(3) 



is proportional to the time variability 6t and the peak volume number density of photons is 
defined as 

n\ = 1 ^^^^ • (4) 



4 TT rsh'^ c Tb Epk 

As shown in Fig. [21 77 pair production within the GRB attenuates the prompt spectrum, 
while electron Compton scattering and 67^6^ have optical depths less than 1 for the 
typical parameter sets we are adopting. The optical depth for 7-ray absorption through pair 
production is the integral above the synchrotron self absorption energy E'^^^ and below the 
Klein Nishima limit e nergy E'f^M = '^e I'emax where the maximum electron Lorentz factor 
is 7 



eKN 



;g^f^. fiDai fc Lull2002h 



Tgrb 



E max2 



dN 



7 jrpl 



IS 



(5) 



In Eq. Othe cross section for 7-7 pair production apair 

opa^r = ^ ^thil - P') (2 /3 i(3' -2) + {3-/3' logiil + /3)/(l - (3))) 

with 



E'y dcmb 



(6) 



(7) 
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in Eq. ([5]) is the number density of photons per energy in the comoving frame arising 
from synchrotron and higher energy processes. The total number density of photons is given 
as a superposition of different synchrotron and possibly higher energy power law spectra 

dN'^ _ dNi'^ dN^'^ dN3\ dN^ 

dE'^ ~ dE'^ ^ dE'^ ^ dE'^ ^ dE'^ ^' 



where 



7 \oi 



for E',,a < E'^ < E' 

pk 



dE'^ E'pk E'pk 

dN^ = ^/J^)/3 f ,<E'< E'ma. 

dE\ E'pk E'pk ^^'^ ^ 

dN^' n' E' 

7 -C/ pfc2 -C/ pk2 



dE'^ E'pk2 E'pk2 



7 ^ -C/ pfc2 • 



(9) 



E'pk and E'pk2 are the synchrotron peak energy and the higher emission peak energy, re- 
spectively. E' asa is the synchrotron self absorption energy, E'max is the synchrotron cut off 
corresponding to the maximum energy the electrons are accelerated to by the Fermi mech- 
anism and E'xN the IC Klein- Nishima limit. E'^in is the IC-boost of the synchrotron self- 
absorption frequency. The ratio R of the high energy component energy flux, {E'pk2Y 

and the synchrotron energy flux, {E'pkY at the peaks is given by 

R=C^, (10) 
Epk 



which is adopted as the most optimistic value i? ~ 10 in our calculations (see above for more 
discussion). Here a and P are free parameters, which we assume to be a = —1 and (3 = —2. 
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Energy (GeV) 

Fig. 2. — Optical depths for processes involving 7-rays in the GRB fireball. Inverse Compton 
and 7 — e pair production do not attenuate the photon spectrum inside the GRB source. 
The only process attenuating the photon spectrum is 7-7 pair production. The attenuation 
is less efficient for higher bulk Lorentz factors and for shorter variability times. The case of 
a GRB having isotropic luminosity 10^^ erg/s are time variability 0.01 second for different 
bulk Lorentz factors is plotted. 
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Prompt photons are absorbed above the cutoff energy Ecutoff 

ml c^Fb^ 

Ecutoff = , (11) 

corresponding to the photon threshold energy for pair production with lower energy photons. 
The prompt photons are not absorbed above the so called thinning energy, the photon energy 
corresponding to optical depth r^.^ = 1 

_ 3ALisoaTHml 

^thinning ~ -poxj- 2 V"*-^/ 

where A = log[2 {2E^ssaE^f'^]/{me Tfe). For high enough bulk Lorentz factors and for short 
enough time variabilities, GRBs are optically thin at all energies. TeV energy photons 
emitted by some processes within the GRB will thus be able to leave the source and be 
observed by experiments like Milagro. 



3. Photon absorption off IR photon fields 



High energy gamma rays are attenuated also when travelling to us because they form 
pairs in collisi o ns wi th low energy photons from the meta-galactic radiation field. Following 
Kneiske et al.l (120041 ). the optical depth of gamma rays depends on the redshift z and the 
gamma ray energy E^ is parametrized by 



Ukg^JE^, Z) 



a. 33 



{EJE, 



\3/2 



(13) 



where Eq = 90 GeV. The attenuation rate versus energy is plotted for different redshifts in 
Fig. [3l The continous lines show the parametrizatio n given i n Eg. [T3l whe r eas th e dashed 
lines represent the best fit model by iKneiske et al.l ( 2004 ). IStecker et al.l ( 20061 ) 's recent 
calculations of intergalactic gamma-ray absorptio n mak ing use of new Spitzer and GALEX 
data diverge from those of iKneiske fc MannheimI (120051 ) at the higher redshifts due to the 
recent discovery that active star formation was taking place in young galaxies at redshifts 
out to beyond 6. This gives larger optical depths at the higher redshifts than previously 
thought. 
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Fig. 3. — Attenuation of 7-ray spectra through 77 absorption at different redshifts. The 
co ntinous hues show th e parametrization of Eq. [13] and the dashed hues the best fit model 
of kneiske eTaP tmi ). 
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4. Prompt and scattered emission 

Throughout the rest of the paper we will assume that both the synchrotron and the 
higher energy GRB spectrum have spectral indices a = —1 and f3 = —2. Also we will 
consider GRB having time variability of 1 seconds and duration of 20 seconds and a bulk 
Lorentz factor of 316. 

Inside and outside the GRB, the spectrum is assumed to be only attenuated by 77 
reactions and the flux which leaves the source and is observed at a luminosity distance 
Dl{z) is 

and Tgrh-yy and Tbkg^-y are the optical depths for pair production inside and outside the source. 



where 
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Fig. 4. — Synchrotron and higher energy prompt flux from one GRB having isotropic lumi- 
nosity 10^^ erg/s, bulk Lorentz factor 316 and time variability 1 second. The slopes a = —1 
and (3 = —2 for both synchrotron and IC spectra. The ratio between the synchrotron peak 
flux and the TeV emission peak flux is assumed to be 10. The solid lines represented the 
fluxes attenuated by 77 absorption inside and outside the sources which are located at red- 
shifts 0.1, 0.3 and 1. The dashed lines indicate the fluxes after attenuation inside the GRB. 
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Outside the GRB, due to 77 interactions with cosmic infra-red background photons (IR) 
and the cosmic microwave background photons (CMB), most of the high energy photons pro- 
duce high-energy pairs. We assume that each lepton of the pair shares half of the energy 
of the initial GRB photon, E^, = The electron-positron flux from the synchrotron and 

the IC fluxes from one GRB at a distance Dl{z) in the observer's frame is 



dA dEp dt 



grb 



E.. 



g — ''"gri)77(2-Be) M g — r(,t;g (2 _Ee ,z)j 



pk 



(16) 



The pair flux will be enriched as long as the burst lasts. The time integrated electron-positron 
fluxes are given by 



dN, 



-dAdEp 



dt' 



dNp 



■'e grb J ^ Ee dt ^^j, 

The pairs will inverse Compton scatter off CMB and IR photons and produce secondary 
photons, which will in turn interact with IR photons and generate other pairs. Multiple 
inverse Compton scatterings will happen until the energy of the secondary photons is no 
longer sufficient to trigger a pair production with the IR photons. When the energy of the 
scattered photons is insufficient to produce a subsequent pair, the photons will travel to 
the detector without undergoing further absorption. When simulating the series of inverse 
Compton scatterings we assume that all interactions happen very close to the source, at the 
redshift of the source itself. We also assume that the magnetic field is stronger than 10~^^ 
Gauss in order for the flux to be isotropically radiated. The number of secondary photons 
per unit volume per photon energy interval created in the vicinity of a GRB through one 
inverse Co mpton scatter ing of the pairs in Eq. (fT6l) off cosmic of the background radiation 
photons is iGaissen (Il990l ) 



(17) 



iliciE^, z, Liso))^. 



dN^ 



3'-^ UE^dV 



I grb 



dEp I de. 



-cmb' 



daic 
dE^ 



i,E^, Ep dcmb) ^cmbi^^cmbi z) 



dNp 



'dA dEp 



igrb 



where '^'^'c{E^^E,ecmh) either the differential Thomson cross secti on (for low energ y 



pairs) or the differential Klein-Nishina formula (for high energy pairs) (ISchlickeiserl |2003| ). 
Ucmbi^cmb, z) is the dens ity of GBR photons per unit volume per photon energy interval at 
redshift z (jPeeblesI Il976l ) . If inverse Compton scattered photons have a high enough energy, 
they will be absorbed by the IR photons. 



In order to evaluate the contribution from scattered 7-ray emission from all GRBs we 
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input the observed GRB redshift distribution (iGuetta et al.ll2005l : iFirmani et al.ll2004j ) 

ncRBiz) = ^ / ^^^^ (^-«) '^o(^..o) . (19) 

In Eq.ffTOl). $o(-^iso) is the luminosity function, defined as the comov ing space densi ty of 
GRBs in the interval log{Liso) and log{Liso) + dlog{Liso)- Following ISchmidtl (119991 ). we 
derive the luminosity function $o(-^iso) by assuming a broken power law with lower and 
upper limit, Liower = ^"^q- and Lupper = lOLjso*, respectively, where Lisa* = 4.410^^ erg/s. 



The normalization constant a 



ISO 

^ Lisa* ' 



for Lio 



upper 



for Liower ^ LigQ <i Lupper 



(20) 



-(1 

71 ^ 



V) + i(-l + ^: 



72 ^ 



(21) 



is obtained by imposing the integral of $o(-^jso) to give unity. We assume 71 = —0.1, 72 
5i = 30 and 62 = 10. 



The integration over Lj^o is performed over the interval indicated in lGuetta et al.l (120051 ). 
In Eq.(fT9l) the redshift distribution of GRBs Rgrb{z) is 



where 



Rgrb{z) 



22 



o3Az 



(22) 



(23) 



(1 + z)3/2 

and PGB.B = 0.4:4:Gpc~^ yr~^ is the observed rate of GRB per differential co-moving volume. 

dh(i + zfda^dn 



dV 

dz 



VLmil + zf + Vtkil + zf + Vti^ 



(24) 



where da is the angular diameter distance and dh is the Hubble distance. 



The secondary source function in Eq.(|T8l). iterated over multiple scatterings, is inte- 
grated over the comoving line of sight distance dl 



dl = dh 



dz 



+ z? + nk(i + z? + nA 



(25) 
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and weighted over the number of GRBs tlgrb per unit time to obtain the total counting rate 
of the diffuse background flux 



scattered - ^^^^^^^ 



^^^"^^'^'-^^ = j dl j dz (g,c(^7.^.^-o)),,,nGi?,B(^)e-^-«-(^-^) (26) 



where the integration over the redshift is performed from redshift 0.01, corresponding to 
about 40 Mpc, to redshift 10. As we mentioned earlier if we assume that the intergalactic 
magnetic field is stronger than 10^^^ Gauss the beamed emission from GRBs is isotropically 
radiated. So the flux contribution for each burst is lower by the beaming solid angle divided 
by 47r. However there are more GRBs, increased by 4 7r divided by the beaming solid angle, 
which contribute to the diffuse background. So the two beaming factor corrections cancel 
out and the unknown beaming fraction does not influence Eq. fl26l) . 



The total flux from the prompt GRB emission is given by the following integral 

dN^ , r dK 

^7 "'^ grb 



p _ Iprompt _ I 1 prompt -, ( ^\ „-nkg^AE^,z) /cy^-. 



where 



dE^dA Jq dE^dtdA 



dt{ ^ yyj (28) 



is the differential gamma ray prompt flux from each GRB in Eq. (IT^ , measured in GeV ^cm ^ , 
^GRBiz) is the redshift dependent number of GRBs per unit time given in Eq. (fT9l) . 



In Fig. [5] we plot the sum of the prompt and scattered GRB emissions assuming an 
average value of the bulk Lorentz factor equal to 316. The average time variability of the 
bursts is assumed to be 1 second and the duration 20 seconds. The input energy flux at the 
high energy peak Epk2 is assumed 10 times stronger than that at the synchrotron peak Ep^, 
i.e. i? = 10 (see Fig. [I]). 

For higher average bulk Lorentz factors Fb, the pair production inside the GRB source 
region is less efficient and therefore the 7-ray flux escaping the GRBs will be higher and 
both the prompt and the scattered fluxes in Fig. will be higher. If the average time 
variability which is related to the dimensions of the GRB shock radius, is lower, then the 
number of photons emitted through electron synchrotron emission inside the GRB source 
region is higher and consequently we expect higher prompt and scattered fluxes. The second 
peak energy Epk2 is the energy at which we expect the maximum TeV energy flux and, as 
pointed out, is constrained by the observed synchrotron spectrum to be between 10^ and 
10^ times the synchrotron peak energy Ep^. If we consider a lower value of Epk2i we have 
more photons at lower energies which scatter off the CBM and IR photons and therefore 
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the cascade processes started by the GRB photons off the interstellar radiation fields will 
be interrupted sooner and will produce lower scattered photon fluxes. A lower second peak 
energy Epk2 might also imply a more efficient pair production absorption inside the GRB 
source region, as shown in Fig. [21 Finally assuming a shorter average time duration T for 
the GRBs less photons will be injected by the GRBs and the scattered flux will be lower. 

As a back of the envelope calculation, by assuming that the average amount of energy 
released by each GRB is of the order of 10^^ erg, the GRB rate is one per day and the 
average redshift is 1, the total flux emitted by all GRBs is roughly of the order of 5 x 
10~^ GeV cm'"^ sr~^. This is roughly consistent with our more rigorous calculations. 



5. Conclusions and discussion 

Current limits for the contribution of blazars and other sources to the extragalactic 
diffuse emission indicate that between 25 and 50 per cent of the extragalactic diffuse emis- 
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from GRBs should contribute to the diffuse extragalactic emission, too. In fact, outside the 
GRB source, due to interactions with cosmic infra-red background photons, most of the high 
energy GRB photons produce produce high-energy electron-positron pairs. The pairs inverse 
Compton scatter off CMB photons and produce electromagnetic cascades. In this way the 
beamed GeV-TeV GRB emission is re-processed and converted into an isotropic MeV energy 
emission, which contributes to the extragalactic diffuse emission. 

We have modeled the emission from GRBs as due a low energy synchrotron spectrum 
and a higher energy IC spectrum. The IC spectrum extends up to Klein-Nishima limit. 
Also we have assumed the higher energy IC emission from GRBs to be 10 times stronger 
than the lower energy standard synchrotron emission. This assumption is equivalent to 
requiring that the GR B inverse Compton Y-paramater is equal to 10. Using the BATSE 
peak flux distribution iGuetta et al.l (120051 ) derived the GRB luminosity function assumed 
for our estimates. The luminosity function can be approximated by a broken power law with 
a break peak luminosity of 4.4 x lO^^ergss"^, a typical jet angle of 0.12 rad, and a local 
GRB rate of 0.44/t^ Gpc"^yr~^. Using the redshift and luminosity distributions derived by 



Guetta et al.l (120051 ) we have summed up the prompt and scattered emissions from all GRBs 
in the universe. The sum of the prompt and scattered emission from all GRBs is shown 
in Fig. O for a particular choice of GRB parameters. From our optimistic model the 7-ray 
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emission from prompt and scattered GRB emissions can provide only a small fraction of the 
extragalactic diffuse emission at EGRET energies and cannot therefore explain the missing 
flux. From Fig. Othe sum of prompt and scattered emission from synchrotron and 10 times 
stronger IC emissions from all GRBs is less than the part of the extragalactic diffuse emission 
not explained by blazars, admitting that the blazars explain between 25 and 50 percent of 
the isotropic extragalactic emission. Our optimistic GRB model does not overproduce the 
extragalactic 7-ray background. We can thus conclude that TeV energy prompt fluxes from 
GRBs can be at least 10 times bigger than the synchrotron flux without violating the limit 
imposed by the extragalactic diffuse emission. In allowing TeV energy fluxes from GRBs of 
compar able or even stronger int ensity to the MeV fluxes, our result is consistent with those 
models (jPe'er fc WaxmanI 120041 ) which predict fluences at TeV energies similar to those at 
MeV energies, the observations of high energy emission from GRB 970417a with M ilagrito 
dAtkins et allboosh and the analysis of GRB 941017 done bv iGonzales et all J2003h . which 
pointed out the existence of a second flux from GRB 941017 that cannot be explained as 
synchrotron emission. Experiments like Milagro and future missions like GLAST or mini- 
HAWC should therefore devote part of their efforts to investigate VHE emissions from GRBs. 



We note that the recent Swift detection of GRB 060218 (ICampana et al.ll2006l ) suggests 
that there rnight be a low-luminosity population of GRBs with a much higher event rate 
( iLiang et al.l 120061). On the other h and, these bursts have low luminosities and tend to be 
X-ray flashes (ICampana et al.ll2006l ). which would compensate their high event rate. Future 
analyses are needed to reveal the contribution from this category of GRBs. 
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Fig. 5. — GRB Total Emission from bursts having average bulk Lorentz factors 316, time 
variability 1 second and duration 20 seconds. The slopes a = —1 and jS = —2 for both 
synchrotron and IC spectra. The ratio between the synchrotron peak flux and the TeV 
emission peak flux is assumed to be 10. The flux is attenuated by 77 absorption inside and 
outside the sources. 
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